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Abstract 

We study the steady state of the two-species Asymmetric Simple Exclusion Process 
(ASEP) with open boundary conditions. The matrix product method works for the deter- 
mination of the stationary probability distribution. Several physical quantities are calculated 
through an explicit representation for the matrix products. By making full use of the relation 
with the continuous big g-Hermite polynomials, we arrive at integral formulae for the parti- 
tion function and the «-point functions. We examine the thermodynamic limit and find three 
phases: the low-density phase, the high-density phase and the maximal current phase. 



1 Introduction 

Nonequilibrium behaviors of physical systems are one of the most important theme of the statis- 
tical physics. Although the general theory is partially absent, a lot of efforts have been made to 
study several kinds of simplified models that possess nonequilibrium properties. One of the most 
studied model is the Asymmetric Simple Exclusion Process (ASEP). This is a model of many 
particles on a one-dimensional system. Due to an exclusion interaction, the model essentially 
belongs to a many-body system, however, the model is solvable and deep analysis can be made 
by means of the matrix product method. There are many variants for the ASEP. The two-species 
version of the ASEP exhibits remarkable behaviors in many ways, such as construction of a 
shock profile [1-3] and spontaneous translation breaking [4-7]. 
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One of the most interesting findings for the ASEP is the relation with the theory of the q- 
orthogonal polynomials [8-12]. In this paper, we extend the relation for the two-species version 
of the ASEP with open boundary conditions. We will see that the continuous big g-Hermite 
polynomials play a central role for this model. By the help of the ^-calculus, we calculate the 
bulk quantities and the rc-point functions for the steady state and obtain integral formulae for 
them. The integral formulae are useful for evaluation in the thermodynamic limit. We clarify the 
existence of the boundary-induced phase transitions. 

This paper is organized as follows. In Sec. [2]the model is introduced and the matrix product 
method is presented for the two-species ASEP with open boundaries. In Sec. [3] a preliminary 
for the g-calculus is given. Several formulae for the continuous big g-Hermite polynomials are 
presented. In Sec. |4]an explicit representation for the matrices and vectors in the matrix product 
method is provided. This representation is the key for understanding the relation between the 
model and the continuous big g-Hermite polynomials. In Sec. |5]an application of the formulae 
in Sec. |3]is made to the calculation of the quantities for the steady state of the model. We obtain 
integral formulae for the partition function and the n-point functions. By taking the thermody- 
namic limit, phase transitions are appeared. A comparison with the result for the one-species 
ASEP is made. The last section is for conclusion. 

2 Model 

The ASEP is a one-dimensional stochastic process of many particles. Take a one-dimensional 
lattice of size L with two boundaries. Consider two species of particles, denoted by 1 and 2 
for the first- and second-class particle, respectively. Let denote an empty site. We assume an 
exclusion interaction of the particles such that there is at most one particle, 1 or 2, on each site. 
Particles hop in a random way on the lattice. The rates of hopping to the right and left nearest 
sites are fixed asymmetric, say p R and p L . By scaling time, we put p R = 1 and p L = q. That is, 
the dynamics of the particles in the bulk is 

10 *± 01, 20 i 02, 21 i 12. (2.1) 

q q q 

We mainly study the case < q < 1 since the case q > 1 is mapped by the parity transformation. 
The symmetric case q = 1 is special. Since we are interested in the nonequilibrium behavior, we 
postpone the symmetric case. At the boundaries, the injection and removal of the particle 2 are 
allowed. At the left (resp. right) boundary, the injection rate is a (resp. 6), and the removal rate 
is y (resp. J3) if the target site is empty. We take a,/3, y, 6 positive. If we choose ft = 8 = for 
example, the system will encounter jamming of particles. We do not consider such an extreme 
situation. We assume that the injection and removal of the particle 1 are not allowed, and hence 
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the number of the particle 1 is fixed to N. The system is illustrated in Figured! If the particle 
1 is absent in this system, we have the one-species ASEP. We remark that the dynamics of the 
particle 2 in the bulk is irrespective of the presence of the particle 1 (see Eq. (12.1D ), while the 
injection and removal of the particle 2 at the boundaries can be affected by the particle 1 due to 
the exclusion interaction between the particles 1 and 2. 

Let Tj = 1,0 and cr, = 1 , be the particle numbers at site i for the particles 2 and 1 , re- 
spectively. The stochastic process is defined as the time-evolution of the probability distribution 
P({r,-, cr,}; t) with respect to the particle configuration {t ( , ct,} = (t 15 • • • , t l , ctj, ■ • ■ , <t l ). We use 
the bra-ket vector form in the basis of l(^)^ =1 |t,, ct ( ))|t; = 1, 0, a, <= 1, o}. The master equation is 
described by 



-\P({T i ,o- i };t))) = H\P({T i ,o- i };t))). 
at 



The evolution operator is defined as 



(2.2) 



L-\ 



H = V H u+ i +H 1+ H L 



i=i 



(2.3) 



with local operators 
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(2.4) 



(2.5) 



We are mainly interested in the steady state of the system. The stationary probability distri- 
bution can be constructed by the matrix product method. We prepare the matrices D, E and the 
vectors (W\, \V) satisfying 
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Figure 1: The two-species ASEP with open boundaries. The black circle denotes the particle 2 
and the white one denotes the particle 1 . Only the particle 2 can get in or get out at the boundaries. 
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(2.7) 



' L 



with some D,E,A. The choice D = I,E = -I, A = O is easy to treat. This gives the relations 



DE - qED =D + E, 
AE-qEA =A, 
DA - qAD = A, 
(W\(aE -jD) = (W\, 
(J3D - 8E)\V) = \V). 

Then, we have the steady state in a matrix product expression, 

L 

P {{tu crt)) = -?-(W\ Y\(TiD + o-iA + (1 - Ti)(l - cn)E)\V), 



(2.8a) 
(2.8b) 
(2.8c) 
(2.8d) 
(2.8e) 



(2.9) 



i=l 



where Z L N is the normalization constant. We often call the normalization as the partition func- 
tion. It is convenient to consider the grand canonical ensemble. Introducing the fugacities £ and 
£ for the particles 2 and 1, respectively, the steady state is generalized to 

l ( 

P({Ti,Oi};££)= — l —-(W\r\{T^D + a^A + {\ - t ; -)(1 -adE)\V), (2.10) 



where the partition function is simply written as 



Z L ({,D = (W\({D + £A+E) L \V). 



(2.11) 



The case £ = gives the one-species ASEP. In fact, the ^ ^ 1 case is achieved just by modifying 
the boundary parameters and the fugacity hence in the below we often consider the £ = 1 case 
and omit £ without any confusion. The canonical partition function is reproduced by an integral 
around the origin, 



(2.12) 
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We can obtain physical quantities for the ASEP as follows. We denote the expectation value 
of with respect to the steady state by (0). For the grand canonical ensemble, the particle 
densities and their variances are calculated as 

Pi = <T t ) = |^logZ L (£<D| f=1 , Ap 2 2 = (|^J logZ^,^, (2.13) 

Pi = (o-i) = ^logZ L (£,a Ap] = li-jL] logZ L (£,0. (2.14) 

Since the particle hopping is set asymmetric, there is a particle current even in the steady state. 
The particle current is calculated by the matrix product method as follows. We have the particle 
current of the particle 2 in the grand canonical ensemble as 

J 2 = Prob{- • • 20 • • • } - <?Prob{- • -02 • • • } + Prob{- • • 21 • • • } - gProb{- • • 12 • • • } 
= ^jz-AmC-'iDE - qED + D£A - q£AD)C L ~'~ l \V) 

z L (D ' 

where we have used (I2.8al) and (I2.8cl) . Similarly, in the canonical ensemble we have 



(2.15) 



h = (2-16) 

We have the particle current of the particle 1 in the grand canonical ensemble as 

7j = Prob{- • • 10 • • • } - <?Prob{- • -01 • • • } - Prob{- • -21 • • • } + gProbj- • • 12 • • • } 

= —{W\C~\^AE - qE£A-D£A + q^AD)^' 1 ' 1 \V) 
Zl 

= 0. (2.17) 
Obviously, we have Ji = in the canonical ensemble. 



3 Continuous big g-Hermite polynomials 

The ASEP is closely related with the theory of the ^-orthogonal polynomial. In the two-species 
model considered here, the continuous big g-Hermite polynomials play the central role. In this 
section, we give a preliminary for the ^-calculus exploited in the analysis. For more information 
and summaries on the ^-calculus and the ^-orthogonal polynomials, please refer to [13-15]. Note 
that formulae appeared in this section hold for \q\ < 1 . 
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We recall the multiple ^-shifted factorial for n = 1, 2, ■ 

s n-l 



(flu-- - ,a s ;q)„ = [ | |~[(1 -a r <?*)- 



(3.1) 



r=\ k=0 

For n = 0, it is defined to be 1. The basic hypergeometric function is defined as 

(a u --- ,a r ;q) k /. .,, nw^-r 



r(ps 
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(q,bu--- ,b s ;q\ 



((-)V*>) 



(3.2) 



In terms of the basic hypergeometric function, the continuous big g-Hermite polynomial is de- 
fined as 



H n (x;C\q) = C n 3<p 2 



q- n ,C^\C^' ie 
0,0 

(q; q) n 



^ (q; q) k (q; q) n - k 



q;q 

(£j e ;q) k e^ e 



(3.3) 
(3.4) 



with x = cos 9. Note that H n (x; £'\q) is a ^-orthogonal polynomial with respect to x of degree n 
as a descendant of the Askey-Wilson polynomial: 



H n (x;C\q) = P n (x;{',0,0,0\q). 
The three-term recurrence relation is known as 

H n+l (x; C\q) + £q n H n (x; C\q) + (1 - q^H^x; £\q) = 2xH n (x; C\q) 
with = and H = 1. The orthogonality relation of //„(.*;; reads 
^dfl(<7,e 2ie ,e- 2ie ;<7) 



Jo 



# m (cos 9; £\q)H n (cos 9; ?\q) = (<?; q) n 5„ 



/0 2^'^,^;^ 
We introduce supplementary two-variable polynomials 

<2„(x,y; t u h) = (ht 2 xy; q\t^\<p 2 
In particular, putting t\ = A and t 2 = in 2„(^c,y; we have 

(<?; q) n 



(3.5) 



(3.6) 



(3.7) 



q n ,hx,hy 




t x t 2 xy,Q 


q,q 



(3.8) 



F n (x,y;A) = ^ 

fe=0 



(^^yV"*, 



(q; q) k (q; q) n . 



(3.9) 



and this polynomial satisfies the recurrence relation 

F n+1 (x,y; A) + Axyq n F n (x,y; A) + (1 - <20xyF„_!(x,y; A) = (x + y)F n (x,y; A) 



(3.10) 



with F_i = and F = 1. Obviously, we have # n (cos 0; = F„(e ie , e~ ie ; f ). 

The bilinear sum formula for Q n (x, y; t\, t 2 ) holds as follows. When t\t 2 xy = s\s 2 zw, 



oo 

y — 1 



«=() 



(<?, s\S2Zw;q) n 



Q„(x,y;t u t 2 )Qn(z,w; Si, s 2 ) 



(jtyxyz, Tt 2 xyz, ts y xzw, ts 2 xzw; q)^ 
(rt\t 2 x 2 yz, txw, Tyz, rxz; q)™ 
x & W 7 [rtit 2 x 2 yz/q; hx, t 2 x, s { z, s 2 z, rxz\q, ryw] , 



(3.11) 



where a very-well-poised r+1 r series is shortly written: 

r+l W r (aua 4 , ■ ■ ■ ,a r+l \q,z) = r +i<f>r 



1/2 1/2 

U\, qciy , — qciy , ci^, • • • , @ r+ \ 



1/2 1/2 , , 
a 1 ' , -a/ ,qai/a 4 ,-- - ,qa\l<x 



r+l 



<?,z 



(3.12) 



If we put x = e 1 ", y = e~ 10 , z = e 1 ^, w = e" 1 ^, Eq. (13.1 II) is reduced to a bilinear sum for 
the Al-Salam-Chihara polynomial. The proof of (|3.1 II) is in the same line as [16]. We need 
specializations of (13.111) . Setting t\ = S\ = £, t 2 = s 2 = and x = e" 9 , y = e~ 1<9 , z = e 1 ^, w = e^ 
gives 



V -#„(cos 6- £\q)H n (cos 0; ?\q) 

^ (<?; qk 



n=Q 



( Te i(e+,*) ; re i(e-« ; Te -i(e+^) 5 TC -i(e-« ; ? ) 



"302 



Ce ie ,£e- ie ,r 



(3.13) 



To obtain a symmetric expression, we have used a transformation formula 

{e/ a, de/bc; q) 



302 



a, Z?, c 


de 




' a£>c 



"302 



J, de/bc 



q, e/ a 



(e, de/abc; q) 

for generic a, b, c, d, e. Next, setting t\ = £ ,t 2 = s { = s 2 = and x = e 10 , y = e" 16 * gives 



(3.14) 



V -f/„(cos 6- £\q)F n (z, w; 0) 

^ (9; <?)« 



ra=0 



(t£z, t£'w; q) c 



(rze ie , rze~ ie , rwe. ie , rwe~ ie ; q) 
Here we have used Heine's transformation formula. 

(az; q)oo 



"202 



^e ie ,re- ie 



q,r zw 



(3.15) 
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a, b 
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q,bz 
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(3.16) 
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With the help of the relation (13.61) . another kind of bilinear sum is obtained: 
V -H n (cos 9; C\q)H n+ i(cos 0; g\q) 



n=0 



2(cos - t cos &)(t, T_Ce\ tC^\ q)o 
J - r 2 [ (jeW+®, TeW-®, Te-^\ re"^"^; q) K 

(jqeP+V, TqeW-®, rqer^ 9 ^, Tqe'W-®; q% 



"3^2 



q,r 



"302 



q,rq 



(3.17) 



In order to take the limit r — > 1 in (|3.17l) . we apply the following transformation formula for the 
first term of RHS of (T3T71) : 



302 



a, b, c 
d, e 



de 



abc 



(e/b,e/c;q) c 
0, e/bc;q) x 



"302 



d/a, b, c 




d,bcq/e 


q,q 



(d/a, b, c, de/bc; q) m 

H 30? 

(d, e, bc/e, del abc; q)^ 



e/b,e/c, del abc 
de/bc, eq/bc 



q,q 



(3.18) 



As a result, we obtain 

H n (cos6;t;'\q)H n+1 (cos(f);£'\q) 



n=0 



(<?; q) n 



= (cos - cos 6) Uq, Ce w , ?<T W ; qU + (q, f e 1 *, C^; q)J /(e i(e+ «, e i( ^ , e^, ^\q\ 



r'e ie ,re- ie ; 



<T' r 



(q,£'qe l0 ,£'qe- w ;q)c 



(6+<f>) np i(9-</>) (7P -i(e+0) ne -i(8 



"302 



q, ce*, 

Cqc ie , £'qe- ie 



q,q 



(3.19) 



4 Representation for the matrix products 

In the analysis of the steady state of the ASEP, full use of the matrix product method can be made 
if we have an explicit representation of the matrices and the vectors satisfying (|2.8I) . In fact the 
matrices are closely related to the g-boson operators. If we set 



D = 



1 



1 - 



-d+d), 



E = 



1 



1 - 



<l+e), 



A = 



1 



l-q 



{de - ed), 



(4.1) 



then d and e are the annihilation and creation operators of g-boson, respectively, and satisfy the 
relation 



de - qed = 1 - q. 



(4.2) 
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On the basis of the g-deformed Fock states, the explicit representations of the matrices are as 
follows: 



d = 



^JT^ 
^fl^ 




e = 












V 1 ^ 

o V 1 -<i 2 



(4.3a) 



A = diag(l,g,g 2 ,---)- 
The boundary vectors are determined from (|2.8dl) and (|2.8el) . yielding 

(W\ = (w ,w l ,---f, |V> = (v ,Vi,---), 

where 

w n = F n (a, c; 0)/ ^(q;q) n , 
v„ = F n (b, d; 0)/ ^(q; q) n 

Here we have introduced more useful boundary parameters: 



where 



K% y = -^[(l-q-x + y)± V(l -q-x + y) 2 + 4xy] . 

For a,fi, y, 6 positive, the parameters satisfy the inequalities a > 0, b > 0, -1 < c < 0. 
0. 

Set = (1 - q)C,. We introduce a vector, called the bulk vector, 

IV(cos 0)> = (/* (cos 6; C), ^i(cos 9; £), h 2 (cos 6; £), ■ ■ ■ f , 
where the nth element is 

fc„(cos 6; C) = H„(cos 6; g\q)j ^(q\q)n- 

From the three-term recurrence relation of the continuous big g-Hermite polynomial 
have the eigenvalue equation 



(4.3b) 

(4.4) 

(4.5a) 
(4.5b) 

(4.6) 

(4.7) 
-1 < d < 



(4.8) 

(4.9) 
(EH), we 



(d + £A + e)\hAx)) = 2x\hAx)), 



(4.10) 



In the terminology of the orthogonal polynomial, the matrix (d+£'A +e) is the Jacobi operator for 
the continuous big g-Hermite polynomial. We define (/?^(cos ff)\ as the transpose of \hg>(cos 6)). 
The orthogonality relation (13.71 ) guarantees the completeness of \hp(x))\ 

r d6 (q, e 2ie , e~ 2ie ; 

1 = X s feg ^ 1 *'^'"^ (4 ' n) 

Equation (|4.11D shows that the spectrum of (d + £'A + e) is {2 cos 0; 9 e [0, 7r]} = [-2, 2]. 

In fact, the representation given here is only valid when \a\, \b\, \c\, \d\ < 1; otherwise, the 
inner products of the boundary vectors and the bulk vectors would become divergent series and 
ill-defined. Therefore, the quantities of the ASEP should be calculated first in the well-defined 
region of the boundary parameters and then analytically continued to other regions. An explicit 
representation available for all parameters is not found yet. 



5 Nonequilibrium phase transitions 

Using the explicit representation in the previous section, we derive integral formulae for several 
quantities of the two-species ASEP. Evaluation of the integrals is done in the thermodynamic 
limit L,N » 1 . 

We begin with the partition function. First, we set £ = 1 and assume \a\, \b\, \c\, \d\< 1 . By the 
use of (14Tb , (14.101) and (14.111) . we have 



Z L (0 

= (W\(2 + d + £A + e) L \V){\ - qY L 
Jo 2n(^ e ,^;q) 



(W\(2 + d + £A + e) L \h?(cos 9))(h ? ,(cos 9)\V)(l - q)~ 



2n (^'e'U'e-^), 



'■(W\h ( , (cos 6)}(2 + 2cos6) L (h^(cos6)\V)(l - q)~ 



(5.1) 



The inner products in the integral are calculated from (I3.15I ). In this way, we arrive at the follow- 
ing integral expression: 



Zl(0 



f 

Jo 



dQ 



(q,Ca,Cb,Cc,Cd,e Il0 ,e- 2l0 ;q) 



2n (£'e ie , £'e- ie , ac ie , ae~ ie , be ie , ber®, ce ie , ce~ ie , de ie , de~ ie ; q) 



x i4>7 



q, ac 



X[2(l + cos#)] L (l -q) 



202 
L 



Z'b,£d 



q, bd 



(5.2) 
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Rewriting this expression with z s.t. z + z 1 = 2 cos 6 yields 

dz (q,Ca,Cb,Cc,Cd,z 2 ,z~ 2 \q)c 



Z L (0 = 



4mz (£'z, £'/z, az, a/z, bz, b/z, cz, c/z, dz, d/z; q) 



X 2 02 



q, ac 



202 



?b,£'d 



q, bd 



xKl+zXl+z- 1 )]^!-^)" 1 . 



(5.3) 



Here the integral is along the unit circle. Now, we consider other regions of the boundary parame- 
ters. Since the integrand has a number of poles that move depending on the boundary parameters, 
analytic continuation of the integral is made by deforming the integral contour such that no pole 
crosses the contour while changing the parameters. Therefore, the contour should always en- 
close the poles eq k and always exclude the poles (eq k )~ l for k e Z> and e = a,b, c, d, '. Let 
M e = maxfm e Z> ; eq m > 1). If we suppose f(z) has no pole at z = eq m , (eq m )~ l , m = 0, • ■ • , M e 
for e = a,b, c, d, then the contour integral can be decomposed into two parts; the sum of the 
residues at z = eq m , (eq m )~ l , m = 0, • • • , M e and an integral along the unit circle, 

dz (z 2 ,z~ 2 ;q)oo 



Aniz (£'z, C lz, az, a/z, bz, b/z, cz, c/z, dz, d/z; q) 

M <> < 2„2m „-2„-2m. 



f(z) 



(a l q lm ,a- z q- £m ;q) c 



(q-' n ; q) m (q, a 2 q m , ^'aq m , £'/aq- m , abq m , b/aq~ m , acq" 1 , c/aq' m , adq m , d/aq~ m ; q\ 
+ [a «-» b, c, d, £'] 

dz (z 2 ,z~ 2 ;q)oo 



-f(aq m ) 



+ 



§■ 



-f(z). 



(5.4) 



/unit circle 4mz (£'z, £'/z, az, a/z, bz, b/z, cz, c/z, dz, d/z; q)oo 
With this decomposition, the integral expression (15.31 ) is valid for all parameter regions. Intro- 
ducing the fugacity for the particle 2 with { 2 requires the replacement a — > {~ l a, b — > {b, c — > 
{~ x c, d — > {d, —> C{~ 1 ■ Th e grand canonical partition function is therefore obtained as 



dz 



(q,£'r l a,Z'b,£T z c,Cd,z\z-*;q) 



4mz (Cr l z, £T x lz, t l az, {~ l a/z, {bz, {b/z, t l cz, fc/z, {dz, {d/z; q\ 



X 202 



C{' l z,C{~ l lz 
£'{- 2 a,C{- 2 c 



q,{ ac 



202 



C{- l z,C{- l /z 



q,{ bd 



x[(i+^z)(i+^- l )] L (i-^r L . 



(5.5) 



We can see that £ = reproduces the one- species partition function in [1 1] up to a factor. 



Let us examine the thermodynamic limit, large L and N = Lp\. It is difficult to evaluate the 
integral (15.51) . but the asymptotic behavior for large L is obtained relatively easily. We should 
make a careful remark. The relation between the fugacity £ and the particle density p { is given 
by (12.141) . Unfortunately, the one-to-one correspondence is not valid for all parameters. Thus, 
we should also switch to the canonical ensemble for some cases. 
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First, we suppose that £' > a,b and £' > 1. Recall that = (1 -q)£. The decomposition (15.41) 
shows that the leading term is from the pole z = £ ', and of order [(1 + £')(1 + g 2 /£')] L . In 
this case, we have p\ = from (12.141) . Note that p t is a positive increasing function of £' and 
satisfies pi — > as £' — > 1 and pi — > 1 as — > oo. Hence, the one-to-one correspondence is valid 
in this parameter region. If we put £o = yz^> m e condition for the one-to-one correspondence is 
equivalent to £ > a, b. The canonical partition function is just obtained by the relation 

z L (e,D^z L , N (fx N . (5.6) 



For other parameter regions, the one-to-one correspondence of pi and f fails since the N = 
part of the grand canonical partition function becomes dominant and the N ^ part does not 
appear in its asymptotic form. Now, suppose that a > £o> o. > b. In the case where the number of 
the particle 1 is fixed to N, the partial sum is expanded as 



2m£ 



= V Res r" N+1) Z L (f, CM + V Res r mi) Z L (f, 0# 



(g Wzg l ,i;b/zq c/zq ,1jd/zq l \q)i{z 1 \q) o 
(q' e , q~ e /z 2 ; q) e (^- i azq' e , £bzq- { , ^~ x czq~ l , %dzq~ l \ q\ 



X 202 



+ [z 



-t -til 

q ,q lz L 

J;~ l alzq~ l ',g~ l c/zq~ 
X[(l+tz)(l+Z/z)] L (l-q) 

l/zl 



-L 



202 



q ,q lz 

gb/zq- £ ,t;c/zq- £ 



q, g 2 bd 



(5.7) 



Integrating with respect to z with the decomposition (15.41) gives the leading term from the pole 
Z = {^~ x aq~ l y x for a small t. Note that the second term of (15.71 ) gives the same value as the 
first term. Due to the factor q m , the leading term of Z LiN is the term taking both I = in the 
summation and the residue at the pole z = (£ _1 a) _1 . The same argument holds for the case where 
b > £o> b > a. 

With these evaluations of the partition functions, the particle current and the particle densities 
are calculated from (|2.13D - (|2.16I ). As a whole, we have the following behavior in each parameter 
region: 



(A) a > fa, a > b 



N-L 



{^a z ;g)M -q) 

(q, tpb/a, c/a, ^ 2 d/a\ q) c 



-a" w [(l +a)(l +fa)] L , 



(5.8) 
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1 + a (1 + a) l L (1 + a) 1 

(B) b > £ ,b > a 



P2 = 



l+b 



-Pu 



(l+b) 2 L' 



j 2 = a-g) 



(l+b) 2 



(5.11) 



(C) a, b < £ 



Z LN (t 2 ) - zgjtf 2 ) 



((„f;?Wi-?) 



Af-L 



rtl+^oXl+f/^, (5-12) 



(a/^o,f-b/^o,c/^o,^ 2 d/^ ;q)c 



1 - Pi t 1 _ P? 1 _ P? 

P2 = — r 1 , Ap^ = — fi, / 2 = (l- 9 )__n. (5.13) 



2 f " 4L 

It is interesting to observe that 



Z £?^ 2 ) - (f> ^TT Z S^)- (5-14) 



With respect to the behavior of the particle 2, we have three phases; just as the one-species 
ASEP with open boundaries, which is recovered by putting pi = in the two-species model. 
Accordingly, the three phases are the phase A: the low-density phase, the phase B: the high- 
density phase and the phase C: the maximal current phase. See Figure [21 By checking the 
discontinuity of the derivatives of the partition function, the phase transition is found to be of the 
second order between the phases A(B) and C and of the first order between the phases A and B. 
On the critical line between the phases A and B, the low-density phase and the high-density phase 
coexist. Indeed, a shock profile of the density of the particle 2 is expected on this coexisting line. 
In the phase A, the density of the particle 2 is not suppressed by the existence of the particle 1. 
On the other hand, in the phase B, the density of the particle 2 is suppressed by pi, and in the 
phase C, suppressed by pi/2. The existence of the particle 1 also causes the suppression of the 
maximal value of the current of the particle 2 from (1 - q)/4 to (1 - g)(l -p 2 )/4 and hence shifts 
the lines of phase transitions between the phases A and C, and between the phases B and C. 
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Figure 2: The phase diagram of the two-species ASEP. There are three phases; the phase A: the 
low-density phase, the phase B: the high-density phase and the phase C: the maximal-current 
phase. The double line indicates the transition of the first order and the solid line indicates the 
transition of the second order. 



The physical explanation of the shift of the critical line between the phases A(B) and C is as 
follows. As mentioned before, the dynamics of the particle 2 is unaffected by the presence of the 
particle 1. However, due to the exclusion interaction between different species, the rates of the 
injection and removal of the particle 2 at the boundaries have effective values, depending on the 
density of the particle 1 near the boundaries. The phase diagram may be the same as the one for 
the one-species, drawn with respect to the effective boundary parameters. Instead, with respect 
to the original parameters, it appears that the critical lines are shifted. 

We can explicitly obtain the n-point function in an integral form. Similarly to the partition 
function, we insert the completeness relation (14.1 II) n + \ times in the matrix products, and arrive 
at an integral formula. Let rjj = r, or cr,. The result is as follows: 

n+l r 



1 



1 



Z L (\-q) L 

*2<p2 



ti-ri p. 

Di 



-[(l+^Xi+z; 1 )]'' 



-if-i-i 



x (az\,alz\,cz\,clzi,bz n +\ 



q,ac 



Amze (Cze, £'lzt\ q% 

CZn+l > £ /Zn+1 

> b/z n +i,dz n +i,d/z„+i',q)J 



202 



q, bd 



x H ^i k =Ti7(z k ,Zk + \) +Xm k =<ri§(Zk,Zk + \)\ , 



(5.15) 



where i ■ 

FALSE 

excludes the poles at ze = {zt-\q k )~ l ,{z~ e \q k )~ l ,{£,' q k )~ x (k 6 Z>o). The contour C\ encloses the 



and i n+ \ = L + 1 in the products, and^ denotes the characterstic function: ^true = 
: 0. The contour C ( (■£ = 2, • • • , n + 1) encloses the poles at ze = Zc-\q k , z^ t q k , Cq k an d 
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poles at zi = aq k , cq k and excludes the poles at z\ = (aq k y 1 ,(cq k y 1 (k e Z> ). The contour C„+i 
also encloses the poles at z n +\ = bq k ,dq k and also excludes the poles at z n +i = (bq k )~ l , (dq k y l 
(k e Z>o). By the use of the summation formula (13.191) . (h(cos + d)\h(cos 9k+i)) is summed 
up as 

J(Zk,Zk+i) =4mzk+i— -r, ; — o{z k+ \ - Zk) 



(<l>zl +v Z k l{,q) 



1 



+ ^fc+i + l lzk+v -z k - l/zk) [(q, Czk, C/Zk, q)oc + (q,Czk+i,C/z k +i;q)o 



x (ZkZk+i,Zk/zk+i,Zk+i/zk, l/zkZk+\',q)' c 
C {q,q£zk,qC lz k ;q) oo 



"302 



q, £zk+\,C lzk+\ 

qCZk, qC'Zk 



q,q 



2 (qZkZk+i, qzk/zk+i, qzk+i/zk, q/zkZk+i ; q) 
By the use of the formula (13.131) with r = q, (/z(cos ^)|^'A|/z(cos 6*^+1)) is calculated as 



(5.16) 



S(zk,Zk+i) 



C(q, qCzk,qC '/z k ;q)c 



(QZkZk+i , qzk/zk+i, qzk+i I Zk, ql ZkZk+i ; q)c 



"302 



q, Czk+\,C /zk+i 

q£'Zk, qCzk 



(5.17) 



It is difficult to evaluate the n-point functions. The phase transitions and the density profiles for 
the two-species ASEP with specialization q = y = 6 = were discussed in [17, 18]. In fact, in 
the phase A and the phase B, phase separation occurs. An intuitive explanation is that if there 
are many particles 2, by the exchange rule 21 ^± l q 12, the particles 1 are pushed from right to 
left and are accumulated near the left boundary. By the particle-hole symmetry, if there are a few 
particles 2, the particles 1 are accumulated near the right boundary. 



6 Conclusion 

We have studied the two-species Asymmetric Simple Exclusion Process (ASEP) with open 
boundaries. For the steady state, we have calculated the partition function through the matrix 
product method and obtained integral formulae. We examine the asymptotic behavior of the 
integral in the thermodynamic limit. Accordingly, we have found three phases depending on 
the boundary parameters, just as the one-species version of the ASEP with open boundaries. We 
have furthermore calculated the n-point functions and obtained integral formulae. In the analysis, 
the ^-calculus is heavily exploited and the continuous big g-Hermite polynomials is important. 
To extend the relation between nonequilibrium particle models and the theory of ^-orthogonal 
polynomials may be a fascinating future problem. 
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